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^ Abstract 

^ A link between structural ordering and slow dynamics has recently attracted much 

attention from the context of the origin of glassy slow dynamics. Candidates for such 
structural order are icosahedral, exotic amorphous, and crystal-like. Each type of 



o 



tj^ order is linked to a different scenario of glass transition. Here we experimentally 

o 

^ access local structural order in polydisperse hard spheres by particle-level confocal 

g microscopy. We identify the key structures as icosahedral and FCC-like order, both 

^ statistically associated with slow particles. However, when approaching the glass 

Q transition, the icosahedral order does not grow in size whereas crystal-like order grows. 

O 

I— I It is the latter that governs the dynamics and is linked to dynamic heterogeneity. This 

questions the direct role of the local icosahedral ordering in glassy slow dynamics and 
suggests that the growing lengthscale of structural order is essential for the slowing 

CN 

down of dynamics and the nonlocal cooperativity in particle motion. 
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Upon cooling, a liquid usually undergoes a first order transition to an ordered ground 
state (crystal or quasi-crystal). However, it is generally possible to avoid the transition and 
to form a metastable supercooled liquid. Upon further cooling the dynamics slows down 
dramatically by many orders of magnitude, leading to the impossibility to equilibrate the 
system in the experimental time scale below a temperature Tg named the glass transition 
temperature. The nature of the glass transition (thermodynamic or kinetic, structural or 
purely dynamical) is still a matter of debate after decades of intensive research (see [H |2] 
for comprehensive reviews). 

A part of the answer is believed to come from the dynamics in itself: a supercooled 
liquid is dynamically heterogeneous (see [2] for a review) and the characteristic size of the 
dynamical heterogeneity is growing when approaching the glass transition [3, 4] (although 
non-monotonic behaviour has been observed recently [5|). The dynamical arrest may then 
be the analogue of the slowing down observed near a critical point when the characteris- 
tic size of the fluctuations diverges, although slowing down at a particle level is absent in 
ordinary critical phenomena. Here we note that the lengthscale defined by the dynami- 
cal heterogeneity is not static (one-time spatial correlation) but dynamic (two-time spatial 
correlation) . 

Modern spin-glass-type theories of the structural glass transition [BHH] introduce a possible 
static length called 'point-to-set' length that makes neither assumption nor description of 
the 'amorphous order' underlying the glass transition. It was demonstrated theoretically [9] 
for lattice models that the growth of a dynamical length scale imposes a lower bound on 
the point-to-set length. However this constraint is expected to become relevant only in a 
low temperature regime inaccessible to present-day simulations |5]. Indeed static lengths 
have been measured uncorrelated to the dynamic one above this regime [TU]. The relation 
between the dynamical and the static length scales is a corner stone of the debate between 
geometrical frustration-based descriptions, facilitation-based descriptions, and the random 
first order transition (RFOT) theory. 

In order to unveil static quantities, it is tempting to link dynamical heterogeneity to struc- 
tural heterogeneity. Stable structures should move less than unstable ones. Widmer-Cooper 
& Harrowell [11] demonstrated that the initial configuration of a simulation, independently 
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of the initial dynamics, has a role in the heterogeneity of the dynamics. Berthier & Jack [12] 
pointed out later that this influence does not exist at the particle level, but on larger scales: 
medium-range particle configuration can be statistically linked to medium-range dynami- 
cal heterogeneity. However these arguments are based on the iso-configurational ensemble 
(starting several simulations with different dynamics but with the same initial configuration), 
a method that has no experimental equivalent. 

Given that a glass shows no long-range periodic order, structures different from a crystal, 
icosahedral [T3HT5] or amorphous order [6], are most often pointed as responsible for the 
dynamical arrest. Since the pioneering work of Frank [16], icosahedral order is the archety- 
pal model of amorphous structures [17] . Locally, an icosahedron maximises the density of 
packing, but its five-fold symmetry is incompatible with long-range periodicity. Geometri- 
cal frustration associated with icosahedral ordering is proposed to be crucial for vitrification 
[T3HT5] . It was also proposed that the percolation of icosahedral order is responsible for glass 
transition [18]. Evidence of icosahedral ordering has indeed been found in very dense hard 
sphere packings [TOl [T9ti22] , in various simulation models [131 UHl I23ti25] and more recently 
in experimental dense metallic liquids [261 127] and glasses [281 129]. 

Geometrical frustration in icosahedral order (or more generally locally favoured structures 
of the liquid) plays a central role in original versions of the spin-glass-type theory [131 E] 
and the frustration-limited domain theory of the glass transition (see [15] for a review). The 
former theory considers the icosahedral ordering under frustration. On the other hand, the 
latter theory takes as a reference state an unfrustrated icosahedral order existing in a curved 
space, even if (and because) this reference state is impossible for the real system to reach. 
If the space were curved, the locally preferred order of the liquid (icosahedral) would form 
a continuous ordered phase upon cooling. In the euclidean space, however, this transition 
is avoided, yet the avoided transition temperature still acts as a critical point giving rise 
to diverging lengthscales (icosahedral domain size and defect size) that would explain the 
dynamical arrest. This theory suggests that the frustration-limited domains are the origin 
of the dynamical heterogeneity. Both theories are still very popular and icosahedral order 
is often cited as the reference amorphous structure particularly in metallic glasses [261 - 
[29] . Yet, to our knowledge, numerical systems where a link between icosahedral order and 
dynamical heterogeneity has been observed [23ti25] . always have a (quasi-)crystalline ground 
state including icosahedra. Thus, icosahedra could be regarded as a structural motif linked to 
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the crystalline ground state in these systems. From this point of view, these results support 
neither the spin-glass nor frustration-limited domain approaches based on frustration in 
icosahedral order, but rather the crystal-based approach (see below). 

An alternative way of explaining the phenomenology of the glass transition is based on 
the structure of the crystal. Crystallisation, even if avoided, may have an influence on the 
supercooled fluid or the glass [301432] and is taken as a reference state. Simulations pHl25| [33] 
and experiments [33] show that a local bond orientational order compatible with the ground 
state (crystal) symmetry does extend to medium range in several model supercooled systems 
at accessible temperatures. In this class of systems, these transient but relatively stable 
structures seem to be correlated with the slow regions because they are low free-energy 
configurations. Furthermore, they lack extended positional order and thus could not be 
detected by macroscopic diffraction experiments. We emphasize that high crystal-like bond 
orientation regions are not crystal nuclei with solid nature, but transient medium-range 
structural order formed in a liquid as thermal fluctuations. 

In the present work, we focus on the colloidal (hard spheres) supercooled liquid [3l] for 
two reasons. First, hard sphere-like colloidal particles can be tracked by confocal microscopy, 
giving access to the positions and dynamics of individual particles [351 136], thus allowing mi- 
croscopic analysis of the dynamical and structural heterogeneity. Second, hard spheres have 
a well defined crystalline ground state of face-centred-cubic (fcc) or hexagonal-close-packed 
(hcp) symmetry and locally favoured structures of icosahedral symmetry. This enables us 
to distinguish the different scenarios of glass transition, in particular, the frustration-limited 
domain scenario, the amorphous order scenario, and the medium-range crystalline order 
scenario. We will show that both kinds of order are statistically associated with slower than 
average particles. However icosahedral order remains local, whereas the crystal-like order 
correlation length grows when approaching the glass transition. Moreover, the medium- 
range crystalline order is associated with the slow regions of the dynamic heterogeneities. 
We confirm that both structural and dynamical length scale grow in the same (critical-like) 
way. Our result suggest that, when not compatible with an avoided crystallisation, local 
ordering seems to play a minor role in the dynamical aspects of the glass transition. 
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RESULTS 



Dynamics 

We were able to follow the volume fraction dependence of the dynamics over nearly 
three orders of magnitude (see Fig. [T]), confirming the basic phenomenology observed in 
experiments [21 IMti36] as well as simulations [33] of similar systems. In Fig. [T^ we fit the 
self-intermediate scattering function by a stretched exponential oc e"'-*/'^"^'' that defines the 
structural relaxation time Tq. Here /3 is the stretching exponent. 

The non-Gaussianity of the dynamics is monitored by the kurtosis of displacements dis- 
tribution (Fig. [T|d) 
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where Ar{t) = \r{to + t) — r{to)\ is the displacement of a particle after a time interval 
t. At a given volume fraction, 0:2 (t) peaks at the characteristic time scale of the dynamic 
heterogeneity f^^lcf)). 

The volume fraction dependence of Tq, is super-Arrhenius and we fit it in Fig. [T]i by 
the Vogel-Fulcher-Tammann (VFT) law = r°exp(D0/ (0o ~ 4>)), where D is the fragility 
index, 0o is the hypothetical ideal glass transition volume fraction, and r° is the microscopic 
timescale. We notice that f^^ corresponds to the end of the plateau, always slightly shorter 
than Ta- 

The spatial heterogeneity of the dynamics is characterised by a four-point susceptibility 
Xiit) and a four-point structure factor 6*4 (g,t) ^7\. To the reader unfamiliar with these 
terms (defined in Methods), we note that 5*4(5, t) ^ J^/^s is the structure factor of a system 
containing only the slow particles. 

As shown in Fig. the small wave-vector behaviour of 5*4(5, t) is well described by the 
asymptotic Ornstein-Zernike function in Fourier space: 

where ^4(t) is the four-point correlation length (see Methods for the details of the fit). We 
found that at each volume fraction the maximum of x^i^) is reached at times longer than 
the relaxation time r^. In summary, we have f^^ < Tq, < t^^. For simplicity, we omit time 
dependence in 5*4, 1^4 and X4 when t = t^^. The volume fraction dependence of the resulting 



correlation length is plotted in the inset of Fig. [Tj In the volume fraction range that we were 
able to study, it follows a power law 

^4 = ao((0O - 0)/0o)-^ (3) 

where = 2/3 and (j)Q is independently determined by the above VFT fit; thus the bare 
correlation length is the only fitting parameter [33]. Note that the lengths are in units 
of the mean diameter a. 

Bond orientational ordering 

The local structures of our suspensions are identified by a detailed analysis of Steinhardt's 
bond orientational order (boo) [13], including original improvements (see Methods). Fig- 
ure |2] summarizes our structure identification results. We find the distinctive signature of 
medium-ranged crystalline order (mrco) of FCC type, without excluding some HCP. We find 
no trace of BCC-like structure. As for the aperiodic structures, icosahedral order stands out 
without any other obvious rivals (dodecahedra are present but can be considered as twisted 
icosahedra) . 

We establish that in our system the tendency towards icosahedral order is best monitored 
by the order parameter called wq, whereas Qq is the natural crystal-like bond orientational 
order parameter (see Methods for the definitions of wq and Qq). The population of local 
structures in the (wq, QG)-'plane is shown in Fig.|2|:;-d, which confirms that the two tendencies 
are present in the supercooled liquid of hard spheres. In particular, icosahedra are present 
even at a relatively low volume fraction (Fig.[2]i). Both tendencies become more pronounced 
for deeper supercooling (Fig. Furthermore, Fig. ^ shows clearly that icosahedral order 
and crystalline order are incompatible and frustrate one another. 

Bond order mobility 

In the introduction we mentioned the dynamic propensity invented by Widmer-Cooper 
& Harrowell [Tlj as a trick only accessible via simulations. In our system, since we know the 
relevant structural order parameters {wq and Qq), we can compute the square displacement 
of all the particles with the same initial local structure (iso-bond order ensemble) and call 
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it bond order mobility. We define the wg-mobility as 

/El|r.(t)-r.(0)f 5«-^6) 
^Aw„t) ^ ( ^ ^ ) , (4) 

i 

where the brackets denote ensemble averaging. At^^Qq, t) and Ar'^{wQ, Qq, t) can be defined 
in the same way. The iso-configurational propensity is very powerful since it does not require 
any a priori information on structures and thus can probe any static order if it exists. This 
feature is absent in the bond order mobility, which instead has the merit of accessibility even 
in a system that cannot be restarted ab libitum from exactly the same overall configuration. 
In general the two quantities are not equivalent. 

Figure [s] shows the bond order mobility (calculated at the characteristic time of the 
dynamical heterogeneity, f^^). Both structural order parameters have a clear correlation 
with the dynamics: the more structured is the environment, the more likely the particle will 
be slow. For all volume fractions, the mobility displays the same linear dependence upon wq 
(Fig. |3]d). Once scaled by the mean square displacement all points collapse on the same line. 
This is not the case for the Qq dependence (Fig. ^) which is more and more pronounced 
with increasing the degree of supercooling. In our most deeply supercooled sample, the 
crystal-like environments can be in average ten times slower than the bulk, whereas even 
the almost perfect icosahedra are in average only 40% slower than the bulk. A complete 
screening of all the possible BOO mobilities excludes a link of slow dynamics to any other 
local symmetries (including exotic amorphous order), at least in our system. 

To explain our mobility data, we can refer to the cage picture, where a particle is consid- 
ered as trapped in the cage formed by its neighbours. Once the particle has escaped from its 
cage, it is free to diffuse. In general, this picture is not correct due to dynamical heterogene- 
ity, or nonlocal cooperativity of motion. In the case of the icosahedral environment, however, 
it seems to hold: the normalisation by the bulk MSD takes into account the diffusion once 
out of the cage, and the remaining universal dependence upon wq holds the information 
about the "quality" of the cage. Efficient packing makes the 13 particles icosahedra more 
stable than a disordered structure. Thus the central particle has a low probability to escape 
from its cage and start diffusing. After cage escaping, however, it diffuses in average like 
any other particle. This suggests that the influence of icosahedral order on the dynamics is 
only local because of its isolated character. By contrast the non-trivial dependence of the 



Qe-Kiobility calls for nonlocal explanations. 
Real space patterns and length scale 

To observe the ordered regions in real space, we define the thresholds Wq and Qg so that 
the bond order mobility at the threshold is half between the bulk and the (extrapolated) 
perfect structure. In our most deeply supercooled sample, this criterion yields Qq = 0.25 for 
crystal-like structures (see Fig. |3]i) and = —0.023 for icosahedral structures (see Fig. [3^). 
These thresholds are arbitrary but coherent with each other. We display in Fig. |4|d typical 
configurations of the ordered neighbourhoods in our deepest supercooled sample. With our 
thresholds, we see only small patches of icosahedral order that are not reaching medium 
range, whereas clusters of crystal-like order are much larger than icosahedral clusters and 
their sizes reach medium range. 

Here it may be worth stressing that MRCO does not correspond to crystal nucleus (com- 
pare Fig.|4^ and b). However sub-critical nuclei are not unrelated to crystal-like BOO. Nuclei 
are systematically embedded into much larger MRCO. One can explain this phenomenon as 
perfect wetting of the nuclei by coherently bond-ordered regions j3B] . We note that the size 
of the crystal nuclei is smaller than the critical nucleus size and thus they appear only tran- 
siently. Furthermore, examples of crystal-like structures without embedded crystal nucleus, 
in addition to the continuity of Qq values (Fig. |2|, suggest that the crystal nuclei are the ex- 
treme part of the bond order fiuctuations. Indeed it has been shown in hard spheres [SHJ-SO] , 
Gaussian-core model [H] and probably Wahnstrom binary Lennard- Jones ^24j that parts of 
MRCO, fiuid regions with a structure reminiscent of the crystal, act as precursors to crys- 
tallisation if the local density there is high enough. However in these systems the presence of 
precursors to crystallisation does not imply the formation of a critical nucleus and growth. 
Furthermore, MRCO is always present (transiently) in the supercooled liquid: MRCO are 
usual fiuctuations of the system, as common as the fast and slow regions of the dynamic 
heterogeneity [SB] . We also note that the average density within MRCO is the same as within 
disordered liquid regions and thus different from within crystal nuclei [38] . 

Now we consider a link between static spatial structures and dynamics. Refiecting the 
short-range nature of icosahedral order, we were not able to extract the associated lengthscale 
that would grow when approaching the glass transition. By contrast, the spatial extent of 
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crystal-like BOO is well described by Gg^r), the spatial correlation function of Qem (see, e.g., 
|33j). Both dynamical and structural lengthscales (defined in Methods) are increasing while 
approaching the glass transition in a coherent manner (the inset of Fig. [T]). We can see that 
near 0o the growth of both dynamic and static correlation length can be well expressed by 
the same power law (Eq. (|3])), although the range is rather limited. In agreement with [33] 
it suggests that the dynamical heterogeneity in hard spheres is a manifestation of critical- 
like fluctuations of the crystal-like bond orientational order parameter, and neither those of 
icosahedral nor amorphous order. In relation to this, we stress that GqIt) is not sensitive to 
aperiodic structures such as icosahedral and amorphous order (see Methods). 

Of course, one can set a more permissive threshold on and observe more icosahedra. 
A percolating network of icosahedral neighbourhood can even be found at a high volume 
fraction if the threshold encompasses enough particles (see Fig. jsj^). However, the imperfect 
icosahedra particles that must be included to form this network have a mobility comparable 
to the bulk (Fig. [3|d). Furthermore, we checked that this percolation is of the 3D random 
percolation class (cluster size distribution follows a power-law of exponent ~ 2.1, see Fig. [3^), 
meaning that a comparable network could be obtained by picking randomly the same number 
of particles (see Fig. |5]i). We conclude that the concept of icosahedral network is not 
physically meaningful at least in our system. 

As a final but an intuitive check, we compare the spatial repartition of stable crystal-like 
order (Fig. ^) and slow particles (Fig. |4]i). As shown in [12], the correlation between struc- 
ture and dynamics is not a one-to-one correspondence at the particle level but a statistical 
relationship on larger length scales. We were able to show visually this correlation by av- 
eraging out short time fluctuations of Qq] by calculating the square displacement of 
each particle over 2Ta ~ A.bf^^ which allow to accumulate many rearrangement events; and 



by coarse-graining this value over the particle's neighbours [12] (see Eq. (14) in Methods). 
Figure |4]i displays the 10% slowest particles according to this criterion and their neighbours. 
The correlation of shape and size is clearly seen between Fig. ^ and Fig. |4ji, although it 
is not perfect. The slight disagreement comes from the following reasons: some clusters 
of FCC-like order disappear soon after the initiation of the displacement measurement and 
some particle groups, which come back to their original positions after the period of 2Ta, 
can be accounted as slow particles. We could flnd no such correlation between icosahedral 
particles (Fig. Eb) and slow regions. 



Note that we recover the correlation between crystal-hke order and slowness by two 
independent methods: directly in instantaneous BOO space via mobility on moderate times; 
and indirectly via positions and shapes of (short time averaged) ordered regions coherent 
with those of slowest regions (on long times). In the former method, quality of statistics 
comes from ensemble averaging over many realisations of each structure. In the latter 
method, it comes from both coarse graining and the length of trajectories that allow to 
accumulate many rearrangement events. 

DISCUSSION 

We have demonstrated that supercooled hard-sphere liquids display local icosahedral or- 
dering as well as structural ordering linked to the avoided crystalline (fcc) ground state. 
Both types of packing have low local free energy and thus are locally favoured. Yet the 
icosahedral order remains local even at deep supercooling and (consistently with the argu- 
ments of P^) is not able to influence significantly the global dynamics. By contrast, the 
slowing down associated with the crystal-like BOO is nonlocal and long-lived. The charac- 
teristic length scale of the crystal-like order grows when approaching the glass transition 
in the same way as that of the dynamical heterogeneity. This suggests that the structural 
origin of the dynamical arrest is linked to the avoided crystallisation and neither to the 
condensation of the local (icosahedral) order of the liquid nor to the exotic amorphous order 
at least in our system. We note that the development of structural order and its link to 
dynamic heterogeneity are at odds with purely kinetic scenarios of glass transition. 

How general is our crystal-like order scenario? Not all glass formers crystallise, indeed 
the crystal structure of many glass formers is unknown. However in a system where the 
crystallisation is possible but unlikely (a definition that include many practical glass form- 
ers), we expect structures locally reminiscent of the symmetry of the crystal to be present as 
common fluctuations of the supercooled liquid. Their relative stability and ability to reach 
medium-range sizes makes them a good candidate to explain the slow regions of the dynamic 
heterogeneity. In that class of systems static and dynamical lengths would grow together. 
Yet the generality and limitation of this scenario need to be checked carefully in the future. 

We note that Charbonneau et al. [10] are not finding growing 6-fold bond order length 
scale in a binary hard sphere mixture. Hopkins et al. |l2t 113] showed that for this size 
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ratio if the two components de-mix they crystalhse independently in FCC or HCP; however 
the symmetry-breaking phase of the mixed system is not known and could possibly be not 
6-fold symmetric. This de-mixing required for crystallization may destroy a link between 
configurations of low local free energy in the liquid state and the crystal structure. Indeed, 
the study of BOO in binary systems with the concept of bond orientational order is not 
promising, as reported in [33| HI] for 2D binary soft disk simulations. 

For some systems, the first order transition avoided upon supercooling leads to a sym- 
metry compatible with icosahedral order. The symmetry-breaking phase can be a quasi- 
crystal [23] or a Frank-Kasper crystal, as recently demonstrated for the Wahnstrom mix- 
ture [2l]. In that case, the icosahedral order can be regarded as part of the defective 
crystalline bond orientational order that slows down the system, so the icosahedral order 
can play a role in the slowing down. However, we showed here that when the icosahedral 
order is not compatible with the crystal symmetry, its role in the slowing down is actually 
minor. 

We do not say that the local (icosahedral) order of the liquid is unimportant in the glass 
transition. Actually it plays a role in preventing the first order transition (crystallisation) to 
happen, thus allowing supercooling and the glass transition |15]. Furthermore, this frustra- 
tion against crystallisation may be linked to the fragility of the glass former [301 132] • Com- 
peting crystal-like and non-crystal-like orderings may be important in many glass-forming 
systems covering colloidal, molecular, polymeric, oxide, chalcogenide, and metallic glasses. 
Indeed, our competing ordering scenario seems very consistent with a recent finding in 
metallic liquids |46j and may explain the high glass-forming ability of bulk metallic glass- 
formers An intimate link between the liquid and crystal structure may also provide a 
mechanism of fast phase change materials [48j . 

METHODS 
Experimental 

We used pmma (poly (methyl methacrylate)) colloids sterically stabilized with methacry- 
loxypropyl terminated pdms (poly (dimethyl siloxane)) and fiuorescently labelled with rho- 
damine isothiocyanate chemically bonded to the pmma. The size distribution, as charac- 
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terised by scanning electron microscopy imaging (see Supplementary Figure SI), showed a 
long tail towards small sizes, leading to the polydispersity larger than 6%. This amount of 
polydispersity allows us to avoid or at least delay crystalhsation [33] (see below the crystal- 
lization monitoring method) but is too low for fractionation to happen |5Q]. The colloids 
were suspended in a solvent mixture of cis-decalin and cyclohexyl-bromide for both optical 
index and density matching. To screen any (weak) electrostatic interactions, we dissolved 
tetrabutylammonium bromide salt, to a concentration of 300 nmol/L |51]. The estimated 
Debye screening length is 13 nm, well below the length scale of the colloids. To avoid errors 
due to the swelling of the particles, their size was measured in the same solvent by the fol- 
lowing technique: short-ranged depletion attraction was induced by adding non-adsorbing 
polystyrene to a dilute suspension to drive particles into contact; the position of the first 
peak of the g{r) gave the diameter a = 3.356 ± 0.03 [im. 

After careful shear melting, the samples are filled into square 500 [im capillaries (Vit- 
rocom). The most deeply supercooled sample showed a slight indication of sedimentation 
even after careful density matching. To avoid this effect we used a thinner 100 [J.m x 1 mm 
cell. The data were collected on a Leica SP5 confocal microscope, using 532 nm laser exci- 
tation. The temperature was controlled on both stage and objective lens, allowing a more 
precise density matching. Our polydisperse colloids were detected using a novel multi-scale 
algorithm (see Supplementary Methods on the details), which allows us to retrieve both 
position and size of each particle with about 1% precision, and this for any arbitrary size 
distribution. We stress that standard particle tracking methods [35l [36] are not able to track 
all the particles in our polydisperse samples (see Supplementary Figure S2). The large size 
of our colloids allows a comparable precision in position for both in-plane and out-of-plane 
coordinates (see Supplementary Figure SI). The particles are tracked in time to extract 
dynamical informations. Typically ^ 5000 trajectories can be followed over a few structural 
relaxation times. 

Dynamics 

To characterise the dynamic heterogeneity, we define a microscopic overlap function 
Wi(t) = Q[b — \ri(t) — rj(0)|], where 9(x) is Heavisides step function, ri(t) is the posi- 
tion of particle z at a time t and b = 0.3(a). The number of overlapping particles is Ns{t) 
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and in the thermodynamic hmit its fluctuations deflne the four-point susceptibihty: 



X4 = i{N^) - {N.)')/{N). (5) 

In a simulation with a flxed number of particles and a limited size, some fluctuations 
are forbidden, thus Eq. ^ do not give the right value of X4 (see [37j). However in our 
experiments we are observing a small portion of the whole colloidal suspension, thus we can 
use time average to sample all the fluctuations of this large system, and Eq. ([s]) becomes 
exact and ensemble- independent. 

The four-point structure factor [37j is 

S,{q,t) = N-\{W{ci,t)W{-ci,t)) - I {W{ci,t)f I), (6) 
where W{q^,t) is the Fourier transform of Wi{t): 

W{ci,t) = ^w,(t)exp(-^q-r,(0)). (7) 

i 

Our experimental data do not have periodic boundary conditions, so we must use a 
window function to ensure the correct correlation, especially at small q. Here we use the 
Hamming window, but we checked that our results are not affected by other reasonable 
choices of the window function. The results are shown in Fig. 

As stressed in [37], the value Si{q = 0) is crucial for the quality of the flt in Eq. (|2|. Here 
we use the value of X4 given by Eq. (|5| and we flt S^^q) for q < 1.5/^4. 

Structure identification 

Steinhardt et al. [T3j deflned the BOO of the i-iold symmetry as a 2£ + 1 vector: 

qim{i) = J7,^ Yi^{e{Yij), (p{Yij)), (8) 
* 

where the are spherical harmonics and Ni is the number of bonds of particle i. In the 
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analysis, one uses the rotational invariants defined as: 



(9) 



21 + 



W£ = ^ 1 j qemiqemiQems, (10) 

mi+m2+m3=0 \ "^1 ^2 "^3 



We =Wi 



\m=—£ 



where the term in brackets in Eq. (10) is the Wigner 3-j symbol. Note that the popular 
E2] We can show a large norm if its denominator (or qi) has small values (particle with low 
£-fold symmetry). Therefore, we have used here we that pulls clearly apart highly symmetric 
particles from the fluid distribution. 

Following Lechner & Dellago [52], we coarse-grain the tensorial BOO over the neighbours: 

Qimii) = ^ yqem{i) + ^g£m(j)j > (12) 

and define coarse-grained invariants Qe and We in the same way as Eqs. ^ and (10). 
Structures with and without spatial extendability are then much easier to tell apart 



For non-extendable local structures like icosahedra, their Qe and We are buried into the 
liquid distribution. With this in mind, we looked for extendable structures in the {Qe, We} 
space, and for non-extendable (aperiodic) structures in the {qe, we} space, with £ = 4, 6, 8, 10. 



In both Eqs. ([8]) and (12), the sum runs over the "neighbours" (labelled with j) of particle 
i. The definition of the neighbours is crucial for the consistency of the analysis. The most 
popular criterion in literature is that particles closer than -Rmax are neighbours, with -Rmax 
the first minimum of g{r). To account for the polydispersity, we used a threshold in the 
scaled distance Vij = rij/{ai + aj) < -Rmax, with -Rmax corresponding to the first minimum 
of g{r). We made a first round of analysis using this definition, constructing maps similar 
to Fig. |2]but noisier. These maps allowed us to discard the possibility of BCC (high Qq, low 
Q4, high Qio) and generally any structure except FCC (high Qq, high Q4 and W4 < 0), HCP 
(high Qq and W4^ > 0) or icosahedron (very negative wq, high q^ and gio). Dodecahedron 
can be found for wq > Wq°'^ ~ —0.00782, i.e., between the liquid and the icosahedral order. 
We are thus left with structures where the particles naturally have 12 neighbours. 
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We then took a new BOO analysis, this time taking the first 12 closest particles (in f) 
as neighbours. This criterion reduces the noise in the distribution of the invariants and 
makes detected structures more easily comparable with crystallographic archetypes. It also 
avoids artefacts near -Rmax in the correlation functions and removes the negative correlation 
between the number of neighbours and the value of the qi. The figures presented here were 
constructed from this second round of analysis. 

Crystal nuclei 

To avoid any confusion between our imperfect crystal-like order and well-formed but small 
transient crystals, we detect crystal nuclei via the most standard method in the literature [191 
[53] . For each neighbouring particle, we compute the normalized scalar product of the (non 
coarse-grained) q^m'- 



The bond between particles i and j is considered crystalline if s{i,j) > 0.7 [55]. A particle 
is crystalline if more than half of its bonds (> 7) are crystalline. 

Polydispersity is an important controlling factor which determines the ease of crystalliza- 
tion of colloidal suspensions (see, e.g., [49J). Polydispersity of our sample was high enough 
to prevent crystallisation in bulk over our experimental time scale. However, heterogeneous 
nucleation from walls sometimes interfered our experiments, thus we had to discard such 
samples. We find no large crystal nucleus in the samples analysed in the text. Even at deep 
supercooling a crystal nucleus never gathers more than 8 particles. This is smaller than 
the critical nucleus size (2-3 particle diameters according to [51] and classical nucleation 
theory). In agreement, crystal-like solid regions appear only transiently without growth, 
i.e. we never observed crystal nucleation in bulk for our samples. Truly crystalline particles 
account for less than 1% of the system. Figure |4] allows to compare the spatial extent of 
crystal nuclei with that of the crystal-like clusters or of the slow regions. It demonstrates 
that trivial crystallisation can be responsible for neither the spatial extent of the dynamical 
heterogeneities nor glassy slow dynamics. 




(13) 
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Slow particles patterns 



The above definition of neiglibours is used consistently in the coarse graining of square 
displacements: 

1 I ^^'^ \ 

^r\{t) = ^^^^ UrMt) + 5^ Ar^(t) 1 , (14) 

where Ar^i(t) = ||rj(t) — rj(0)||^ and the neighbours are defined at initial time. 



Icosahedral percolation 

We studied clusters formed by icosahedral particles as a bidirectional graph. The nodes 
of the graph are all the particles, and at the beginning we have no edge. When a particle 
is 'activated' the bonds between this particle and its 12 neighbours are added to the graph. 
It is often meaningful to activate the particles as a function of an order parameter, for 
example sorted by increasing wq, the most icosahedral particle connects first and we activate 
particles that are less and less icosahedral until percolation and beyond. For our most 
supercooled sample, the percolation threshold is identified near wq ~ —0.011 (see Fig. ^) 
which corresponds to 7.5% of activated particles. We stress that this level of icosahedral 
order corresponds to an average mobility very close (^ 90%) to the bulk (see Fig. [sj^). 

In this way, the cluster size distribution and fractal dimension can be studied system- 
atically at each value of the threshold on the order parameter. In Fig. [5^-b we report 
characteristics of icosahedral clusters when the threshold is set just below percolation (all 
particles with wq < —0.012 are activated). The cluster size distribution shows an exponent 
of 2.1, characteristic of the random percolation universality class. The fractal dimension is 
close to 2. 

One can also activate particles in a random fashion and measure the fraction of activated 
particles needed for percolation to occurs. We performed 100 try of this Metropolis algorithm 
on 160 configurations, in order to estimate the percolation probability function of the fraction 
of activated particles. Results are reported in Fig. |5ji and indicate that the peculiar case of 
icosahedral order clusters has nothing special. 
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Bond order spatial correlation 



We defined the coarse-grained bond order correlation function 



13 E.., ^(n, - r) 



(15) 



whicli is not sensible to aperiodic structures. After baseline subtraction, the envelope of 
GQ{r) is fitted by the Ornstein-Zernike form to extract the crystal-like bond orientational 
order correlation length ^g- We checked that g^ir)^ defined from the non-coarse-grained 



bond order in the same way as Eq. (15), yield similar results. Insensibility against aperiodic 



structures comes from the scalar product in Eq. ( 15 ) [T8] , not from coarse-graining. 
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FIG. 1. Dynamics of the system, a, Decay of the self-intermediate scattering (ISF) function 
computed at wave vector q = 2Tr/ (a). Lines are stretched exponential fits from which we ex- 
tract the structural relaxation time Tq. b, Non-Gaussian parameter, whose maximum defines the 
characteristic time of the dynamic heterogeneity f^^. c, Collapse of the small-g regimes of the 
four-point structure factor on the Orstein-Zernike function (solid curve), d, 0-dependence of Ta 
(circles) and of f^^ (squares). The solid curve is the VFT fit for Tq, {(f>o = 0.600, D = 0.328). Inset: 
Dynamical (^4, squares) and structural (^g, circles) correlation lengths scaled by their respective 
bare correlation length, respectively ^4,0 = 0.206(T and ^6,0 = 0.129cr. The curve is the power-law 
given in the text. All times are scaled by the Brownian time tb- 
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FIG. 2. Population of local particle configurations specified by the boo parameters. 

The value of boo parameter for a perfect structure is indicated by its name's position on each 
map. a in the {Q4,QQ)-plane, b {w4,QQ)-plajie and c-d ('u;6, <56)-pla'ne. a-c are the populations 
at deep supercooling {4> = 0.575 it 0.03); d is the same as c but for a liquid near the freezing 
point {(f) = 0.497 it 0.03). The colour at a given couple of BOO parameters indicate how often 
(log scale) this couple is detected (per units of area in this plane). The arrows stress the ordering 
tendencies: the tendency towards FCC is always visible, a weak tendency towards HOP can also be 
distinguished in b, and the tendency towards icosahedral order (iCO) is visible in c. Conflicting 
tendencies towards FCC and ICO in c is a manifestation of competing orderings, or frustration, 
between them. Straight lines arc drawn at important thresholds: Qq, Wq, Wq°'^ and = that 
differentiates between FCC and HOP. 
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FIG. 3. Bond order mobility, a, Normalised mobility in the {wQ,QQ)-plane for our most deeply 
supercooled sample {(p = 0.575 it 0.03). The colour scale is saturated at 1.5 times the bulk mean 
square displacement, b-c, Normalised mobility for icosahedral and crystalline order parameters 
respectively at volume fraction 0.535 (squares), 0.555 (triangles) and 0.575 (diamonds), all ±0.03. 
Bulk mean square displacement is scaled to be at 1 (horizontal line). Perfect structures are on the 
edge of each plot. The lines are a guide for the eye, stressing the collapse of the i(76-mobility at all 
volume fractions in b and the absence of such collapse in c. The scattering at low volume fractions 
is due to poor averaging of rare structures. Straight lines in a-c corresponds to the important 
thresholds: Q^, Wq and Wq°'^. Examples of crystal- like cluster and distorted icosahedron at the 
respective threshold values are shown in d and e, respectively. 
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FIG. 4. Computer reconstruction from confocal microscopy coordinates in our most 
deeply supercooled sample. Depth is ~ 12a. Only particles of interest and their neighbours are 
displayed. Each particle is plotted with its real radius, a, Particles with more than 7 crystalline 
bonds, b, A typical configuration of bond ordered particles. Icosahedral particles are shown in 
the same colour if they belong to the same cluster. If a particle is neighbouring both crystal- 
like and icosahedral structures, it is displayed as icosahedral. c, Crystal- like particles alone (the 
order parameter was averaged over d, Slow particles with respect to the coarse-grained 

displacement. Due to particles going in and out of the field of view, assignment of particles located 
very near the edges of c and d were not accurate and have been removed. 



25 




FIG. 5. Imperfect icosahedral clusters, a, Cluster size distribution (compatible with random 
percolation) for nearly percolating threshold {wq < —0.012). b, Radius of gyration (indicating a 
fractal dimension df ^ 2 given by the straight line) for the case of a. c, Size of the largest non- 
percolating cluster as a function of the threshold on wq. Percolation takes place -0.011. 
Vertical lines indicate the position of and Wq°'^ respectively, d, Translation of c in terms of the 
fraction of activated particles (dots) together with the probability of the onset of percolation when 
particles are randomly activated (gray shading) . We can see that activating particles as a function 
of icosahedral order and those activated randomly have almost the same percolation threshold 
(Ri 7.5% of the particles activated). 
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